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I. INTRODUCTION 

Quantum graphs and quantum networks have become a popular tool to model quantum dynamics 
in mesoscopic systems in the diffusive regime (Shapiro 1982, Chalker and Coddington 1988, see also 
Dittrich 1996 and Janssen 1998 for recent review articles), as well as, to study spectral statistics 
(Kottos and Smilansky 1997, 1999, Pakonski et al 2001, Ketzmerick et al 2000, Tanner 2001). In 
this paper, the autocorrelation function of the spectral determinant 

N 

Z{d) = e-^(^^+^) det (1 - e'^U) = e-^(^«+^) ^ a„e'"^ (1) 

n=0 

is studied after averaging over a suitable ensemble of quantum graphs preserving the underlying 
'classical' dynamics on the graph. Here, U is a unitary N x N matrix which describes the quantum 
evolution on the graph and the a„'s denote the secular coefficients of the characteristic polynomial. 
The phase factor in front of the determinant ensures that Z is a real function for real 9 after setting 
exp(i(/3) = det{—U) and using the basic property for the secular coefficients 

a„ = e"^a^_„ (2) 

which follows from the unitarity of U. The autocorrelation function of (^) has been studied by 
Kus et al (1993)and Haake et al (1996) after averaging over circular-orthogonal, -unitary and - 
symplectic ensembles of unitary matrices. Related results for ensembles of Hermitian matrices have 
been given be Kettemann et al (1997). All these approaches have in common that the average is 
taken over a set of matrices too large to link the ensemble in the semiclassical limit to a specific 
classical system. In this sense, the results are of random matrix type. There is strong evidence, 
however, that the statistical properties of spectra related to a fixed chaotic classical system already 
follow random matrix theory (RMT). This suggests that the size of the ensembles considered can 
be reduced drastically without loosing universality in the spectral statistics. 

I will approach this problem by defining unitary matrix ensembles corresponding to a specific 
quantum graph. The autocorrelation function of the spectral determinant is studied for these 
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ensembles. It will be shown that the variance of the secular coefficients an follow indeed the 
random matrix result for chaotic classical dynamics up to a m ultiplicative factor. This factor is 
system dependent and will be explicitly derived in section IV This improves the formula given 
by Kettemann et al (1997) by carefully including repetitions of periodic orbits leading to non- 
negligible contributions. Related results have been reported by Keating and co-workers (Keating 
et al 1996) using periodic orbit approximations of the spectral determinant. 



II. QUANTUM GRAPHS AND UNITARY STOCHASTIC ENSEMBLES 

In the following I will use a definition of quantum graphs which is slightly more general than the 
one introduced by Kottos and Smilansky (1997). Consider a graph consisting of N directed edges 
or bonds connecting an unspecified number of vertices. Free, one dimensional wave propagation 
occurs along edges i of length Li. Transitions from an edge i to an edge j is possible at a vertex 
if and only if i is an incoming edge and j an outgoing edge at the same vertex, see Fig. |^. The 
scattering process at a given vertex is described by a scattering amplitude Sy, where i, j label 
again incoming and outgoing edges, respectively. To keep the set of possible networks as general 
as possible, the local scattering amplitudes will not be specified any further other than to fulfill 
basic properties like probability conservation. (One may think of the scattering event taking place 
at a given vertex to be due to interaction with a complex short-range potential or a many-particle 
compound at that point). 




leading to 

coming from 

other vertices 



FIG. 1. Section of a typical directed graph; waves travel along directed edges i,j,k,l,... undergoing 
scattering at the vertices. 

The eigenstates (or resonances) of such a network are given by stationary states of the system. 
These are obtained by considering the N - dimensional wave- vector 4* out where the i-th component 
of ^'out represents the wave amplitude on the i-th edge when leaving the vertex. This vector is 
transformed into a set of incoming waves at the other end of the edge by 

*in = D(fc) *out with dij = e'''^^6ij 

where k represents the wave number. The non-trivial part of the dynamics is contained in the global 
scattering matrix S with local scattering amplitudes describing transitions between incoming 
and outgoing edges i and j, that is, 

*out = S *i„ = S B{k) *o,t = U(fc) (3) 

which in turn defines the propagator U(fc) of the discrete wave dynamics on the network. In 
the following, I will assume that the graph represents a closed system without edges leading to 
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infinity or absorbing vertices. This implies that S and thus U are unitary. Stationary solutions 
are obtained for wave numbers k fulfilling the secular equation 

det(l -U(fc)) = 0. (4) 

Certain restrictions on the topology of the graph and its dynamical properties, that is, on the pos- 
sible transition probabilities |sij follow from the unitarity condition. This point will be discussed 
later in more detail. I will furthermore assume here that the scattering matrix S is independent 
of the wavelength k. 

Statistical properties of quantum systems are linked to the complexity of an underlying classical 
dynamics. For quantum systems acting on a continuous coordinate space with continuous or 
discrete time variable, the classical dynamics revealed in the semiclassical limit is generated by a 
set of Hamilton's equations of motion or for discrete times by a symplectic map. For a quantum 
graph, the wave dynamics takes place on discretised time and space variables and a connection to 
an underlying classical motion is less obvious. The best one can achieve in these circumstances is 
to identify the quantum propagation on a given graph with a Markov process. This probabilistic 
'classical' dynamics is generated by a stochastic transition matrix T with matrix elements given 
as (Kottos and Smilansky 1997, Barra and Gaspard 2001, 2002) 

tij = \uijf = \sijf (5) 

where tij denotes the transition probability to jump from an edge i to an edge j at a given vertex. 
Transition matrices which are derived from a unitary matrix via condition (|^) are called unitary 
stochastic (Berman and Plemmons 1994). A unitary stochastic matrix is also doubly stochastic, 
that is, the rows and columns of T sum up to one. This implies, that T has a largest eigenvalue 
equal to one with corresponding left and right eigenvector ;^(1, 1, 1, ... , 1). Note, that the converse 
is not true, that is, not every doubly stochastic matrix is also unitary stochastic which leads to 
surprisingly complex restrictions on the possible Markov processes which can be 'quantised' in 
terms of eqn. (||), see Pakonski and Zyczkowski (2002). 

What constitutes the semiclassical limit of a quantum graph? The short wave length limit k ^ oo 
for fixed network size N corresponds to a simple scaling transformation, that is, the statistical 
properties of the spectrum are independent of k. Instead, one may construct a 'semiclassical 
limit' for quantum graphs by appropriately increasing the size of the graph. Such a limit should 
be taken in a way which ensures convergence of the stochastic motion towards the dynamics of 
a deterministic map within a given phase space resolution. This semiclassical limit constitutes 
therefore at the same time a classical (deterministic) limit of the stochastic dynamics. Sufficient 
conditions for constructing such a limit will be formulated at a later stage. 

The wavenumber k defines a one-parameter family of unitary matrices U(/c) for a fixed quantum 
graph. A stationary solution exists for /c-values, for which condition (^) holds, that is, whenever an 
eigenvalue exp(i0i(/c)) of U(fc) crosses the real line at -1-1. Correlations in the /c-spectrum are thus 
related to correlations in the eigenphases spectrum (f>i{k) of U(/c) for fixed k. Instead of studying 
the k spectrum directly, one can therefore look at the statistical properties of the eigenphases of 
U(fc) for fixed k and then average over k. 

In the generic case, when the set of lengths Li are all rationally independent, the one-parameter 
family \J{k) covers uniformly a subspace of the unitary group U{N) with the topology of an (at 
most) N dimensional torus. That is, the k average can be replaced by an ensemble average over 
the ensemble Ut of unitary matrices defined as 

Ut = {V e U{N)\V = S D with dij = e''^'S,j, ipi e [0, 2n) Vi = 1, . . . A^} 

where the scattering matrix S is kept fixed and the integration measure for the ensemble is given 

by 

i=l 
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Note that the ensemble average corresponds to averaging over the spectrum of a specific quantum 
system with fixed classical dynamics. All the unitary matrices forming the ensemble are indeed 
linked to the same stochastic process on the graph described by a unitary stochastic transition 
matrix T. I will therefore denote these ensembles as Unitary Stochastic Ensembles IAt (USE) 
(Tanner 2001). 

It has been argued in Tanner (2001) that correlations in the eigenphases spectrum after ensemble 
averaging depend crucially on the spectral gap of T, that is, on the distance of the second largest 
eigenvalue of T from the unit circle. More precisely, it was conjectured that the spectral correlations 
follow random matrix theory in the classical limit if the gap decreases slower than 1/N for N 
oo. I will in the following study a different statistical measure in more detail, namely the auto- 
correlation function for spectral determinants. 



III. THE AUTOCORRELATION FUNCTION FOR THE SPECTRAL DETERMINANT 

The averaged autocorrelation function C{x) of the spectral determinant (Q) may be defined in 
terms of the generating function for the square moduli of the secular coefficients, that is, 

N 

P(z) = ^< |a„p>^,^z", (6) 

n=0 

where the average is taken over a unitary stochastic ensemble as defined above. One obtains 

r(,A < C dez{e + ^)z*{e-^) >u. , . , F(e'") 

< /p d9\Zi9)\^ >Ut ^^^) 

It is in general convenient to rescale the argument such that the correlation function is periodic 
with period one, that is, one considers C{x) — C{^j^). The symmetry property (||) ensures 

< |a„P >=< |aAr_„p > and thus C{x) ^ C{1 - x) . (8) 

Averages will always be taken over USE's unless stated otherwise and I will drop the suffix 
< . >Ut from now on. One obtains immediately 

< a„ >= Vn = 0, . . . , iV 

when averaging over the phases (pi. The variance of the secular coefficients < |a„p > has been 
calculated in the context of RMT by Haake et al (1996). One obtains in particular after averaging 
over the whole unitary group U {N) (CUE), over the ensemble of symmetric unitary matrices (COE) 
or over the ensemble of diagonal unitary matrices (Poisson) 



r 1 CUE 



< |a„P >= < 



^ + N+1 (Q) 

N \ ' ^ ' 

Poisson 

n ' 



Note, that the CUE and COE distributions as well as the CSE distribution not shown here (Haake 
et al 1996) converge to a limiting distribution for large N with respect to the parameter r — n/N 
whereas no such asymptotic limit exists in the Poisson case. It will be shown in the following that 



^It should be noted that a USE is in general not uniquely defined by a unitary stochastic matrix T for A*' > 
2; there are in general several unitary matrices S not connected by multiplication with diagonal matrices 
which corresponds to the same transition matrix T. The statistical properties of all these ensembles are, 
however, expected to be identical in the limit N ^ oo. 
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the ensemble average for USE's without symmetries may deviate from the CUE resuh above. 
Using Laplace expansion of the determinant, one writes 

N 

det (1 - zU) = 1 + ^ z" {-!)" det(U[a]) 

n=l aer{n,N) 

where T{n, N) denotes the set of integer vectors with 

r{n,N) = {ae IN"|1 < ai < as < . . . < a„ < iV} 

and U[a] is the n x n matrix obtained from U by taking only the rows and columns of U with 
indices ai, . . . , a„, that is, 



The variance of the distribution of secular coefficients can thus be written as a double sum over 
all possible subgraphs spanned by the edges a — (ai, . . . , a„), 

< |a„|2 > = < det(U[a])det(Ul'[a']) > 

a,a'er(n,N} 

= Y < det(U[a])dct(Ut[a]) > . (10) 

aer{n,N) 

In the last step, the fact that contributions with a ^ a' vanish identically when taking the ensemble 
average has been exploited. Using Laplace expansion again one writes each of the sub-determinants 
in (10) as a sum over permutations tt of the row (or column) indices, that is. 



det(UH) - ^(-i)"«H. = 5](-i)"n"H.,. 



(t) ■ 

The sum over permutations may be interpreted as a sum over closed (periodic) paths on the 
subgraph a which visit each edge a^, i = 1, . . . n exactly once. These periodic paths will be called 
covering orbits in what follows; they are not necessarily connected. The product of determinants 
in (10) now takes on the form 

a£r(n,N) tt,tt' 

Periodic orbit pairs with tt = tt' will obviously survive the ensemble average and may form an 
important contribution to the double sum. Splitting these diagonal contributions from the off- 
diagonal contributions and assuming that there are no systematic periodic orbit degeneracies due 
to symmetries, one obtains 

1 

= E Per(TH)+ < E i-ir^^'u[aU[a]:, > . 

Here, tla]-^ — |m[q;]^P denotes the product of transition probabilities along the periodic path. 
Furthermore, T[q;] is the sub-matrix obtained from the rows and columns ai , . . . Q!„ of the stochastic 
transition matrix T and per(T[Q;]) is the permanent of this matrix. Note, that the permanent of 
a square matrix A is defined as 



per(A) = EII"*^^ 



(i) 
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A detailed account of the properties of permanents may be found in Mine (1978), Berman and 
Plcmnions (1994) or Brualdi and Ryser (1991). 
One finally obtains for the generating function (^ 

JV 

P{z) — ^ < |a„|^ > = per (1 + zT) + non-diagonal contributions . (11) 

n=0 

The autocorrelation function of the spectral determinant can thus be linked to the permanent of 
sub-matrices of the stochastic transition matrix T. This is in analogy to the connection between the 
spectral correlation function and the traces of powers of T (Tanner 2001). Note, that additional 
contributions to the diagonal approximation have to be taken into account in the presence of 
symmetries as for example time reversal symmetry; I will not consider this case here any further. 

Two questions arise immediately: firstly, what is the range of validity of the diagonal approxi- 
mation, and secondly, how can one extract at least asymptotic results in the large N limit from 
the expression ( |ll|) by for example making a connection to the spectrum of the transition ma- 
trix T? Especially the last point turns out to be quite tricky; dealing with permanents is a very 
hard problem in general (unlike calculating determinants), the effort of computing a permanent 
increases indeed faster than exponential with the matrix size. ^ I will give answers to both the 
points raised above in the next section. 



IV. ASYMPTOTIC RESULTS FOR CHAOTIC DYNAMICS ON GRAPHS 



When studying the autocorrelation function in the limit of large network size, one first has to 
specify, how to take this limit, that is, how to choose families of graphs which may be said to have 
a well defined classical and thus semiclassical limit for — > oo. Taking a very conservative point of 
view, I will demand that an increase of the size of the network has to be done in a way which leaves 
the probabilistic dynamics of the Markov process and thus the periodic orbit structure invariant 
up to a certain cutoff time nc{N); the cutoff time itself needs to increases with N , one can typically 
achieve Uc ~ logA^. Such a definition ensures automatically that typical measures of chaos like 
the topological entropy ht, the K - entropy or the spectral gap converge for large N. Only graph 
families without symmetries and with a finite gap and such exponential decay of correlations will 
be considered in what follows. 

The transition matrices of such families of graphs with chaotic classical limit become increasingly 
sparse for large N; the number of non-zero matrix elements increases typically not faster than N. 
This in turn implies that most of the subgraphs a of size n entering the sum ( p^ ) contain no 
or at least one covering orbit for small n. (I will make no distinction between connected and 
non-connected orbits here.) The diagonal approximation is obviously exact in these cases. Non- 
diagonal contributions will become important for n values for which two or more covering orbits 
exist in a typical subgraph. A transition between these two regimes occurs roughly at times which 
coincide with the radial size of the network, that is, with the mean time to reach each edge from 
every other edge. This transition time nt is of the order 

nt{N) — 

ht 

where ht is the topological entropy. The diagonal approximation can thus be estimated to be valid 
for coefficients < |a„p > with n < nt- 

A connection between the permanent appearing in (|l^) and the spectrum of the transition matrix 
can be established using similar arguments. Consider first the identity 



■^Calculating the permanent of a general matrix is an NP complete problem in the language of algorithmic 
complexity, that is, there is no (known) polynomial algorithm to compute the permanent; see Brualdi and 
Ryser (1991), as well as Mine (1978) for a description of the best known algorithms. 
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det(l - zT) = - z^'Hp) 



(12) 



for determinants, where the product runs over all primitive periodic orbits of the graph, that is, 
connected periodic paths not including repetitions. Here, Up denotes the length of the orbit and 
the weight tp corresponds to the product of transition probabilities along the orbit. One might 
now be tempted to just replace the minus sign on the right hand side in (|l^ by a plus sign, that 
is, to consider the function 



oo 

E 



Z Cn 



(13) 



to obtain a periodic orbit expression for per (1 + zT). The product in ([l^ ) does, however, include 
self-intersecting periodic paths on the graph, that is, periodic orbits which visit the same edge 
more than once. It is one of the magic properties of the determinant, that, after expanding the 
product in (|2|), contributions of intersecting orbits cancel exactly, leaving only terms containing 
orbits and products of orbits, which have no edge in common. The same is obviously not true 
for the expression One can argue, however, that short periodic orbits on the graph are in 

general irreducible, ie, they do not self-intersect (Bogomolny 1992, Tanner and Wintgen 1992). 
Furthermore, taking the classical limit in the sense described above corresponds to enhancing 
the phase space resolution, that is, self-intersecting orbits turn eventually into irreducible orbits 
when increasing N. Again, self-intersections are likely to occur only for times n larger than the 
transition time rif. The coefficients of a polynomial expansion of the permanent in ( [Tl| ) can thus be 
approximated by the coefficients of C+i^) exactly in the regime, where the diagonal approximation 
holds, that is. 



< a„\ 



1 d" 



1 d 



per (1 + zT) 1,^0 « + ■^""^f) 1^=° < 



log TV 



(14) 



The zeta function (|l^ ) can be linked to the transition matrix by the following calculation, ie. 



C+ = + z^'Hp) = exp 



p r—l 



EjfE(-"''*pr-E(-"''^p)") 



exp 

det (1 - z2t(i) 
det (1 - zT) 



(15) 



The matrix T'-^-' is defined by squaring the matrix-elements of T, that is, = tfj. Equation jl^ ) 
is the main result of this paper. It differs from the relation given by Kettemann et al (1997) by 
having a determinant also in the numerator; this term arises by including repetitions of periodic 
orbits in a consistent manner and must not be neglected. It leads to significant contributions in the 



expansion of C+ and thus to the coefficients < 



> in general. These additional contributions 



have also been noted by Keating and co-workers (1996). Their results, derived from semiclassical 
periodic orbit expansions of the spectral determinant, are more general in the sense that they 
apply for generic quantum systems. The resulting periodic orbit formula equivalent to ( [l5| ) could, 
however, not be linked to classical operators. The periodic orbit formulas are notoriously difficult 
to calculate and specific values for the secular coefficients could indeed not be obtained so far. 
Here, one may proceed as follows; after using Newton's method (Kettemann et al 1997) to 



expand the two determinants in (y_5|) in terms of traces, one recovers the coefficients c„ of C,+ in 
eq. ( p^ ) recursively by applying standard formulas for fractions of power series (Gradshteyn and 
Ryzhik 1994). The c„'s converge exponentially fast to a finite value for large n for T matrices 
with a finite spectral gap. The asymptotic behaviour is in particular determined by the residium 
of C-i-(^) at the leading singularity at z = 1. It follows from the relation (|lj) which is valid over a 
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range of n values increasing with N that the coefficients < |arp > converge to a constant value a 
for fixed r n/N and N oo given as 

det (l - z2t(i)) 

a= lim < la,-!^ lim Res^^i — --^ 7^^- (16) 

N^oo ' N~.oc dct(l-zT) ^ ^ 

This result differs from the CUE result insofar that the constant a is in general not equal to 
one. It depends furthermore on the asymptotic limit of the eigenvalue spectra of both T and T^^^ 
and thus on the classical limit of the family of graphs considered. Let me demonstrate this point 
with the help of a simple example. 

Example. Consider the sawtooth map on the unit interval with k legs of slope k. The map has a 
simple Markov partition whose Z-th refinement divides the unit interval into fc' subintervals of equal 
length. The Markov partition defines a transition matrix Tk{N) of dimension N — k^ with k'"'^^ 
nonzero matrix elements all equal to l/k. The so defined transition matrices are unitary stochastic 
corresponding to a graph where each vertex has k incoming and outgoing edges. The spectrum of 
Tfe and T^^^ consist of a single non-zero eigenvalue being 1 and l/k, respectively, independent of 
the order of refinement I. Transition matrices for fixed k form a family with well defined classical 
limit. Following eqs. (|l5|), (|l6|) one expects for the coefficients of the autocorrelation function 

< laop >= 1, < laip >= 1, a =< |a„P 1 - i for 1 < n < ^ 

k 2 

in the limit N ^ 00. (The function C+ does of course not obey the symmetry conditions (|l5|).) 

In the next section, I will consider a non-trivial example in more detail, namely so-called binary 
graphs (Tanner 2000, 2001). 



V. NUMERICAL RESULTS 



The model systems considered are graphs with an even number of edges and transition matrices 
of the form 

^ ( f + <52.+i,) for < z < f , ^ 0, . . . , TV - 1 . (17) 

The graphs represent for dimensions N ~ 2'- the Markov partition of a two-leg sawtooth map with 
complete binary symbolic dynamics. This is a special case of the example above for k = 2. I 
will refer to graphs of the form ([ij) as binary graphs. The transition matrices (|l^) are unitary 
stochastic; one can furthermore show, that binary graphs with dimensions N = p2^ for fixed p 
being an odd integer form a family with a well defined classical limit. The spectral gap is A = 1/2 
for p ^ \ and the spectrum of T (as well as T*^^^) depends on p, but not on the order I. 

In Fig. Ha, the secular coefficients < |a„p > obtained from ( |l0|) are shown together with the 
coefficients of per(l -I- zT) and C+(^) for a member of the p — ?> family with dimension N = 24. 
The coefficients of the various approximations all coincide for n < 5 = log 24/ log 2. The 

coefficients of the permanent function decrease rapidly thereafter, whereas those of the periodic 
orbit product approach a constant (which is 0.625 for this particular family). The coefficients of 
the autocorrelation function still fluctuate wildly for such a small N value. Fig. ^ shows various 
members of the same family for larger TV values. The fluctuations decrease considerably and con- 
vergence to the limiting value a = 0.625 is observed. The asymptotic value depends on the family 
and thus on p; one observes in particular convergence to the asymptotic value a = 1/2 for p = 1 
as predicted in the example in the previous section. 
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FIG. 2. Coefficients of the autocorrelation function for the binary graph family with p = 3; 

a) coefficients < \an\^ > compared with the diagonal approximation obtained from the coefficients of 
per (1 + zT) and the periodic orbit formula C,+ = rip(l + z^^'tp) for A'^ = 24; the asymptotic result for this 
family is 0.625; 

b) convergence to the asymptotic result for members of the same family with larger network size; the 
random matrix result is < |a„p >= 1. 
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